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I NTRODUCT ION 

Adaptive  arrays  based  on  the  LMS  algorithm  [1]  can,  in  principle, 
receive  desired  signals  from  any  arrival  angle.  However,  the  output 
s i gn  a i -  to- i n terf  er ence- p 1 us-no i se  ratio  (SINR)  obtained  from  the  array, 
as  the  desired  and  interference  signal  arrival  angles  and  polarizations 
vary,  depends  critically  on  the  element  patterns  and  spacings  used  in 
t!i  j  array  [21. 

The  purpose  of  this  report  is  to  discuss  a  procedure  that 
may  be  used  to  choose  element  patterns  and  spacings  in  an  adaptive  array. 
The  method  described  is  appropriate  when  the  desired  signal  arrival 
angle  will  not  be  known  in  advance,  so  the  array  must  provide  coverage 
over  some  sector  of  space. 

The  method  is  the  following.  One  starts  with  an  array  consisting 
of  an  initial  set  of  elements,  chosen  by  some  means.  One  evaluates 
the  performance  of  this  array,  to  find  particular  signal  arrival  angles 
and  polarizations  where  the  SINR  is  poor.  One  then  augments  the  array, 
by  adding  extra  elements  one  at  a  time,  to  improve  the  SINR  for  thes  ' 
angles  and  polarizations. 
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This  technique  will  work  properly  as  long  as  the  new  elements 
added  to  the  array  do  not  affect  the  patterns  of  the  original  elements. 
That  is,  the  mutual  coupling  between  the  new  elements  and  the  original 
elements  must  be  small. 

Such  a  method  can  be  used  because,  if  the  original  element  patterns 
do  not  change,  the  SINR  obtained  from  an  adaptive  array  cannot  decrease 
when  additional  elements  are  added.  The  SINR  is  a  nondecreasing  function 
of  the  number  of  elements,  for  all  signal  arrival  angles  and  polari¬ 
zations.  Thus,  once  a  set  of  elements  has  been  found  that  yields  ade¬ 
quate  performance  over  part  of  a  desired  sector  of  space,  one  does  not 
have  to  worry  that  adding  new  elements  will  reduce  the  performance  in 
this  region.  The  new  elements  can  be  chosen  to  improve  the  SINR  in 
other  regions  where  the  original  array  did  not  perform  well. 

The  fact  that  the  SINR  obtained  from  an  adaptive  array  cannot 
decrease  when  more  elements  are  added  is  easy  to  see  on  physical  grounds. 
It  is  known  that  LMS  weights  in  an  array  yield  the  maximum  SINR  that 
may  be  obtained  from  a  given  set  of  elements  for  any  given  set  of  signals. 
(Widrow,  et  al  [1]  have  shown  that  LMS  weights  yield  minimum  mean-square 
error  signal,  and  Baird  and  Zahm  [3]  have  shown  that  minimum  mean-square 


i  <i 


error  weights  are  equivalent  to  maximum  SINR  weights.*)  Suppose  we 
have  an  N-l  element  array,  to  which  we  add  another  element  to  obtain 
an  N  element  array.  Assume  a  given  set  of  signals  is  incident  on  the 
array.  Then  the  SINR  with  the  N  element  array  must  be  at  least  as  large 
as  that  with  the  N-l  element  array,  because  the  N  element  array  can 
always  achieve  the  SINR  of  the  N-l  element  array  by  setting  the  weight 
on  the,.new  element  to  zero.  The  validity  of  this  statement  depends, 
of  course,  on  the  assumption  that  the  patterns  of  the  original  N-l  elements 
did  not  change. 

Since  the  SINR  cannot  decrease  when  a  new  element  is  added,  it 
can  only  increase.  The  amount  of  the  increase,  however,  depends  on 

the  pattern  and  location  of  the  new  element.  In  this  report, 
we  determine  what  condition  the  new  element  should  satisfy  to  maximize 
the  SINR  increase  for  the  case  where  a  desired  signal  and  one  interference 
signal  are  incident  on  the  array.  We  also  give  an  example  to  illustrate 
the  procedure. 


♦Baird  and  Zahm  [3]  prove  this  equivalence  under  the  assumption  the 
signals  are  narrowband.  We  shall  also  make  that  assumption  below. 
However,  in  studying  element  pattern  effects  in  adaptive  arrays,  the 
narrowband  assumption  is  automatically  satisfied,  because  pattern  effects 
are  most  easily  studied  with  CW  signals  anyway. 
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FORMULATION 

Consider  an  N-l  element  adaptive  array.  Let  the  jth  element  of 
the  array  be  located  at  vector  position  7j  with  respect  to  the  origin 
of  a  spherical  coordinate  system  r,e,$.  Consider  first  the  desired 
signal.  Suppose  the  desired  signal  has  polarization  state  Pd*  and 

A 

propagates  into  the  array  from  angle  ecj»4>d-  Let  kd  be  a  unit  vector 
parallel  to  the  direction  of  propagation  of  the  signal.  This  signal 

1.L 

will  produce  a  voltage  x.(t)  on  the  j  array  element: 

J 


V> 


(1) 


where  sd(t)  is  the  desired  signal  waveform  (measured  at  the  coordinate 
origin,  and  expressed  in  analytic  form),  c  is  the  velocity  of  propa¬ 
gation,  and  f j(ed’^d,Pd^  ’s  the  Pattern  response  of  the  jth  element 
to  a  signal  from  angle  ed,<j>d  with  polarization  Pd.  We  may  write  sd(t) 
in  the  form 


d(t)  =  ud(t)e 


j(o)0t+4id) 


(2) 


where  ud(t)  is  the  complex  envelope,  u>0  the  carrier  frequency,  and 
4>d  the  carrier  phase  of  sd(t).  We  assume  \pd  is  a  uniformly  distributed 
random  variable  on  (0,2n). 


*We  shall  not  need  a  complete  description  of  signal  polarization  here, 
but  we  may  think  of  Pd  as  the  point  on  the  Poincare  sphere  [4]  repre¬ 
senting  the  desired  signal  polarization. 
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The  N-l  element  signals  x.{t)  may  be  written  as  components  of  an 

J 

*/*  j'  ^  J 

N-l  dimensional  desired  signal  vector  =  (x-|(t),  X2(t),  --,  xN_i ( t) ) 
(T  denotes  the  transpose).  Under  the  condition  that  the  signals  are 
sufficiently  narrowband,  or  equivalently,  that  the  dimensions  of  the 
array  are  small  enough,  we  have 


u 


d 


#  ud(t) 


(3) 


for  every  element  of  the  array.  In  this  case  the  desired  signal  vector 
may  be  written 


Xd  =  lVt)  e 


j(w0t+*d) 


d  ’ 


(4) 


where  Ud  is  the  vector 


f  1  (0d’^d’Pd^e,- 


J*, 


VfN-l^0d’^d’Pd^e 


Vl  : 


(5) 


/  r-*k.\  th 

and  <f>d  is  the  phase  of  sJt  +  ■  Jc  -  -  J  at  the  j  element,  measured 

J  ^  j*  * 

with  respect  to  the  phase  of  sd(t)  at  the  coordinate  origin: 


O) 

-  -r(?rkfi) 


(6) 
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Next  consider  the  interference  signal.  Assume  this  signal  propa¬ 
gates  into  the  array  from  angle  0.,$.  wi*h  polarization  P^.  The  inter¬ 
ference  signal  vector  will  then  be 


Xi 


j(u>  t+ty.) 

^(t)  e  0  1  Ui  . 


with 


Ui 


^2^e1**i»pi^  G  1 


j4> . 

VfN-llei'*i'Pi>  6  "''I 


and 


ill.  /v 

♦ij  =  (rj*ki)  • 


(7) 


(8) 


(9) 


where  the  notation  is  analogous  to  that  for  the  desired  signal.  u^(t) 
is  the  complex  envelope  and  4>.  the  carrier  phase  of  the  interference, 
and  is  a  unit  vector  pointing  in  the  direction  of  propagation.  We 
have  again  assumed  a  narrowband  signal,  so 


A 


*  «i(t) 


(10) 


for  every  element  j.  Also,  we  assume  4k  is  uniformly  distributed  on 
(0,2u)  and  independent  of  4>d. 


Jim 


.  II 
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Finally,  suppose  that  the  signal  from  each  array  element  contains 

v/'  O 

a  zero-mean  thermal  noise  voltage  n.(t)  of  power  a  ,  and  that  these 

J 

noise  voltages  are  statistically  independent  of  each  other: 


=  0  5«,m 


(11) 


where  E{-}  denotes  the  expectation  and  is  the  Kronecker  delta. 
Also,  we  assume  the  nj(t)  are  independent  of  ^  and  • 

Under  these  conditions,  it  has  been  shown  [2]  that  the  output 
signal-to-interference-plus-noise  ratio  (SINR)  from  an  adaptive  array 
may  be  written 


SINR(N-I)  =  Ed 


uV  -  _ 

dd 


02) 


where  *  denotes  complex  conjugate,  and  £d  and  are  the  input  desired 
and  interference  signal-to-noi se  ratios,  defined  by 


and 


=  7e{ 


u,(tir 


(13) 


5,  •  \ 
a 


(14) 


In  Equation  (12),  we  use  the  notation  SINR(N-l)  to  indicate  that  this 
is  the  SINR  for  N-l  elements. 
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Now  suppose  another  element  is  added  to  the  array,  so  there  are 

N  elements.  (We  assume  no  change  in  the  original  N-l  element  patterns.) 

The  SINR  with  N  elements  will  have  the  same  form  as  Equation  (12),  except 

that  the  vectors  Ud  and  Ui  will  have  N  components  instead  of  N-l.  To 

distinguish  between  the  N-l  element  array  and  the  N  element  array,  we 

let  Ud  and  remain  defined  as  in  Equations  (5)  and  (8),  and  define 
N  N 

new  vectors  and  for  the  N  element  array: 


U"  = 


f  i  ^ed’<*>d*Pd^ 
f2^9d,<{,d,Pd^ 


VfN(0d’<t,d’Pd) 


and 


(15) 


(16) 


N  N 

However,  since  the  first  N-1  components  of  1)^  and  are  just  Ud  and 

N  N  N 

U^,  we  may  write  Ud  and  U-  as  partitioned  vectors.  For  U^,  we  have 


i  . 


(ft  denotes  an  N-1  dimensional  vector  of  zeros)  and 


j*d 

UdN  =  fN^0d*^d*^d^e 


(20) 


It  will  also  be  convenient  to  write  Udf|j  in  terms  of  its  magnitude  and 
phase. 


Jnr 


pd  e 


(21) 


Note  that  nd  may  differ  from  4>d  if  ^(^’^d  ronta1n?'  a  Phase  shift 


Similarly,  for  U.  we  have 


=  I  + 


\oJ  \u 


-  Vi  +  ViN  ’ 


where 


Jn< 

UiN  =  i*4>-j*P-j)e  =  P,e 


Writing  the  SINR  as  in  Equation  (12),  but  with  vectors  and  instead 
of  and  U^,  and  using  the  orthogonality  relations 


T*  T*  T*  T* 

VdN  *  VdViN  *  VdN  =  ViN  =  °’ 


we  find 


SINR(N)  = 


,  t  *  ? 

5d  udud  *  °d 


luX*wJn|2 

'  - 1  T  ★  T 

Ki  +Uilh+pf 
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where 


n  = 


nd  ‘ 


To  simplify  the  algebra  below,  let  us  also  define 


and 

BeJY  =  uju*  . 

Then  Equation  (27)  may  be  written 


(28) 


(29) 

(30) 


SINR(N)  =  Kd 


T  *  ? 

UdUd  +  pd  ‘ 


Be 


JY 


Vi 


.e 


jn 


2  1 


A  +  p 


(31) 


Let  us  first  show  that  SINR(N)  cannot  be  less  than  SINR(N-l), 
regardless  of  Pd,  Pi  or  n  (i.e.,  regardless  of  the  signal  arrival  angles 
and  polarizations  and  regardless  of  the  pattern  and  location  of  the 
Nth  element).  Substituting  Equations  (29)  and  (30)  into  Equation  (12), 
we  have 


SINR(N-l)  =  5d[uT/d  -  . 

Defining 

A  =  [siNR(N)  -  SINR(N-l)]  , 
we  find  front  Equations  (31)  and  (32)  that 


(32) 


(33) 


£ 
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A  = 


Pd  A(A+p2)  +  B2(A+p^)  -  A|BeJ  +  p(jpi  e^n| 
A(A  +  p2) 


(34) 


The  numerator  of  this  expression  may  be  rearranged  into  the  form 


(Apd  ~  Bp . )2  +  2ABPdPn.  p  -  cos(y-n)J  ,  (35) 

which  is  clearly  nonnegative  because  cos(y-n)  <  1.  Also,  since  the 
denominator  of  A  is  positive  by  definition,  we  have  A  >  0  or 
SINR(N)  >  SINR(N-l)  as  stated. 

Now  consider  the  problem  of  choosing  the  pattern  and  location 
of  the  element  to  maximize  SINR(N).  Since  SINR(N)  in  Equation  (31) 

is  unbounded  as  a  function  of  Pd,  let  us  suppose  Prf  is  fixed  and  determine 
the  values  of  and  o  that  maximize  SINR(N).  Maximum  SINR(N)  will 
be  obtained  when  the  quantity 


|BeJY  +  pdpi  eJn| 
A  +  p2 


(36) 


is  minimized  —  i.e.,  when  BeJ  and  p^p^ eJ n  have  opposite  signs.  Hence 
we  first  choose 

n  =  y  -  *  •  (37) 

When  q  has  this  value. 


SINR(N)  =  r 


T  *  O  (B-PhPJ‘ 

♦  °d  -  -T-T-1- 


A+P, 


(38) 


Then  the  value  of  P^  giving  maximum  SINR(N)  is 


With  these  values  for  n  and  p^.  we  find 


SINR(N)  =  Kd 


T  * 

UiUd 


=  5, 


r  ».t 
u 


T  *  T 
N1  N 

d  ud  J  » 


(40) 


which  we  recognize  as  the  maximum  attainable  SINR  if  the  response  of 
the  new  element  to  the  desired  signal  is  Pd.  This  is  the  SINR  that 
would  be  obtained  with  N  elements  when  no  interference  is  present. 

Recalling  the  definitions  of  Pd,  Pi  and  n  in  Equations  (21),  (25) 
and  (28),  we  see  that  the  pattern  of  the  Ntb  element  should  be  chosen 
to  satisfy  the  conditions 


lfN<Wd>  ■  lu>*l  • 


(41) 


and 


J*, 


<fN<WPd>e  -  ”  • 


(42) 


to  obtain  maximum  SINR  for  a  given  set  of  signal  parameters  ©d,  <t>d, 

Pd’  9i’  ^i  and  Pi’  Note  both  the  Phase  ofr  ^(Q»d>,P)  and  the  element 
location  can  be  varied  to  satisfy  Equation  (42). 

Equations  (41)  and  (42)  do  not,  of  course,  define  the  pattern 
of  the  new  element  completely.  They  merely  specify  a  relationship  be¬ 
tween  its  amplitude  and  phase  response  for  two  signal  directions  and 
polarizations.  There  is  considerable  flexibility  left  in  the  pattern 
of  the  new  element  even  after  Equations  (41)  and  (42)  are  satisfied. 

In  a  typical  design  problem,  the  original  N-l  element  array  may  yield 
poor  SINR  for  several  sets  of  arrival  angles  and  polarizations.  One 
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may  be  able  to  choose  the  pattern  and  location  of  the  new  element  to 
satisfy  Equations  (41)  and  (42)  for  more  than  one  set  of  signals  at 
once.  However,  even  if  the  new  element  has  been  optimized  for  only 
one  set  of  signals,  it  will  usually  improve  the  SINR  for  other  signal 
arrival  angles  and  polarizations  as  well. 


AN  EXAMPLE 

Now  consider  a  simple  example  to  illustrate  the  use  of  Equations 
(41)  and  (42).  Suppose  we  have  an  array  of  two  dipoles  located  and 
oriented  as  shown  in  Figure  1.  To  simplify  the  problem,  let  us  consider 
only  linearly  polarized  signals  whose  electric  field  and  direction  of 
arrival  lie  in  the  plane  of  the  paper.  Let  the  angle  e  be  defined  as 
shown  in  Figure  1.  Assume  the  patterns  of  the  two  dipoles  are 

f1(e)  *  cos(0-6O°)  ,  (43) 

and 

f  2 ( 0 )  =  cos (0+60°)  .  (44) 

If  we  define  the  coordinate  origin  to  be  at  the  center  of  dipole  1, 
then 

*d,  =  *i,  *  (45) 

Also,  the  dipoles  are  two  wavelengths  apart,  so  we  have 

<j>d  =  -  4irsin0d,  (46) 

and 

<t».  =  -  4irsin0. 

1 2  1 


(47) 


Using  these  quantities  in  Ud  and  U^,  we  may  compute  the  S I  NR  for  any 
given  ©d  and  0.  from  Equation  (12). 

It  is  shown  in  Reference  2  (in  particular,  see  Figure  9  and  the 
accompanying  discussion)  that  this  pair  of  dipoles  has  a  grating  null 
problem  when  0d=22. 3°  and  0.=39°.  The  lower  curve  in  Figure  2  shows 
the  SI NR  as  a  function  of  0^  computed  from  Equation  (12)  when  0d=22.3° 
(and  for  £d=0  dB  and  ^^=40  dB).  The  effect  of  the  aratinq  null  may  re 
seen.  Interference  at  0=39°  not  only  causes  a  null  at  39°,  it  causes 
another  null  at  22.3°  (the  grating  null).  When  the  desired  signal  is 
near  the  grating  null,  the  SINR  is  low. 

Suppose  now  we  add  another  dipole  to  this  array  to  improve  the 

performance  for  the  set  of  signals  8d=22.3°  and  0^=39°.  Let  the  new 

dipole  have  its  beam  maximum  at  0=0  : 

o 

f3(e)  =  co s(e-e0),  (48) 

and  be  located  a  distance  d  to  the  right  of  dipole  1  in  Figure  1.  The 
problem  is  to  determine  the  proper  values  of  0Q  and  d.  By  computing 
Ud  and  U .  for  the  original  2-dipole  array  in  Figure  1,  we  find  that 
when  0d=22.3°  and  0^=39°,  UdU-=.760  (real).  Hence,  from  Equation  (41), 

0Q  must  be  chosen  so 

|f3(22.30)f3(39°)|  =  |cos(22.3-0o)cos(39-0o)|  =  0.760.  (49) 

Solving  this  for  qq  gives  two  possible  values,  0o=2.8°  and  58.5°. 

Since  the  SINR  achieved  with  the  new  element  included  will  be 
?d  +  Pj  (see  Equation  (40)),  and  since  pd  =  cos(22. 3°-2.8°)  is 
larger  than  pd  =  cos(22.3°-58.5°) ,  we  choose  0O=2.8°  to  obtain  the  largest 
improvement. 
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Next,  Equation  (42)  may  be  used  to  find  d.  Since  the  phase  angle 
of  ^(8)  with  0O=2.8°  is  zero  at  both  22.3°  and  39°,  Equation  (42) 
reduces  to 

^  [sin(22- 3°)  -  sin(390)]  =  -  tt  ,  (50) 

which  yields  d=2A.  Adding  the  third  dipole  gives  us  the  array  shown 
in  Figure  3.  Finally,  computing  the  SINR  for  this  3-dipole  array,  again 
with  9^=22. 3°,  we  obtain  the  top  curve  in  Figure  2.  It  may  be  seen  how 
the  performance  has  been  improved  around  0^39°. 
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